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Abstract 

(*) 

In our previous study, we calculated the transition from factors oi J/ip ^ D^^(^ using the QCD 
sum rules. Based on the factorization approximation, the obtained form factors can be applied 
to evaluate the weak non-leptonic decay rates of J/tp Dj^^^ + M, where M stands for a light 
pseudoscalar or vector meson. We predict that the branching ratio for inclusive non-leptonic two- 
body weak decays of J/ip which are realized via the spectator mechanism, can be as large as 
1.3 X 10^^, in particular, the branching ratio of J/ip D*^ + can reach 5.3 x 10^^. Such values 
will be marginally accessed by the ability of BESIII which will begin running very soon. 
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I. INTRODUCTION 



The decays of J /ip are dominated by strong and electromagnetic interactions via cc an- 
nihilating into intermediate gluons and photon at s-channel. By contrast, the weak decays, 
due to smallness of the strength of weak interaction, are rare processes. Under the spectator 
approximation, one of the charm quark or anti-charm quark in J / ip decays into light quarks, 
and the decay rate of a charm quark (anti-quark) is proportional to G'jpml where Gp is 
the Fermi coupling constant. Numerically the total branching ratio of weak decays was esti- 
mated to at the order of 10~^ Recently, due to remarkable improvements of experimental 
instruments ^nd techniques people turn their interests onto these rare processes from both 

5|, la] sides. The forthcoming upgraded BESIII will be able to 



experiment and theory 



accumulate more than 10^° J/ip per year which makes it possible to marginally measure 
such weak decays in near future. More important, such rare processes are also particularly 
interesting from the viewpoint of theory. On the one hand, it may provide further accurate 
examination of the mechanism which is responsible for the hadronic transition and fully 
governed by non-perturbative QCD effects. One can also expect that such decays may offer 
a unique opportunity to probe new physics beyond the standard model j^, , including the 
minimal supersymmetric standard model, the extra dimension model, the two-Higgs doublet, 
top color-assisted technicolor model etc , in the weak decay of vector mesons. The reason is 
that in such rare decays, weak coupling is rather weak and new physics may have a chance 
to show up. 

In a previous study, we presented a detailed analysis of the semi-leptonic decays of J/ip 
6j, where the branching ratios for such channels were estimated to be at order of 10^^° and 
hence is almost impossible to be observed at BESlll. 

The fundamental ingredients involved in the semi-leptonic processes are the transition 
form factors of J/ib which are evaluated in terms of the three-point QCD sum rules 

(QCDSR) 10, lUl, in that work. Obviously, even though while deriving the form factors 
our goal was to estimate the branching ratios of semi-leptonic decays, under the factorization 
approximation, they can be applied to study the non-leptonic decays. Thus, we will take 
a step forward to investigate the exclusive non-leptonic decays with focusing on two-body 
processes. 



2 



In this study, we will explore the non-leptonic decays J ^ D^*J^ + M, where the final 
states contain a single charmed meson and a light meson M, such as vr, K, p, K* etc.. These 
weak decays are realized via the spectator mechanism that one of the charm quark (anti- 
charm quark) acts as a spectator. In the Standard Model, at the quark level, the Feynman 
diagrams for charm quark decay are depicted in Fig. [H The anti-charm quark decay can be 
obtained analogously by exchanging c ^ c. The effective theory for hadronic weak decays 
have been well formulated jisj]. The most difficult work is to calculate the hadronic matrix 
elements which are governed by the non-perturbative QCD dynamics. 

The non-relativistic QCD can simplify the picture by phenomenologically handling some 
non-perturbative QCD effects and has been widely applied to study some decay modes where 
heavy quarkonium are involved. However, it does not help much for the heavy-light mesons 
where relativistic effects may be significant. 

The first order approximation for the derivation is the factorization hypothesis, where 
the hadronic matrix element is factorized into a product of two matrix elements of single 



currents 
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191]. In this scheme, one element can be written in terms of the 



decay constant of the concerned meson while the other is expressed by a few form factors 
according to the Lorentz structure of the current and meson (a pseudoscalar or vector, for 
example). The non-factorizable effects are incorporated into the effective coefficients which 
are usually assumed to be universal and determined by experiment (Only in some cases, they 
are pertubartively calculable. In reality, these coefficients depend on the concrete processes 
and differ case by case, but the variation may be not very drastic). For the weak decays of 
heavy mesons, such factorization approach is verified to work very well for the color-allowed 
sub-processes. It is reasonable to believe that this conjecture would be valid for J/tp, at 
least for the processes where the color-allowed sub-processes dominate. Thus, the study on 
two-body non-leptonic decays offer an ideal ground to testify the factorization hypothesis in 
the heavy quarkonium system and this test may be more appealing than in decays of D^*^ 
because J/ip contains two heavy constituents. Moreover, they are of great importance to 
discriminate various theoretical tools for the evaluations of transition form factors. 

The structure of this paper is as follows. In section II, the factorization approach for 
the non-leptonic decays is introduced and the formulations are given. In section III, after 
displaying the inputs involved in this work explicitly, the numbers of branching fractions for 
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FIG. 1: Quark diagrams for non-leptonic weak decays of J /"ip. (a) represents the color allowed pro- 
cesses; (b) represents the color suppressed processes; (c), (d) represent single-Cabibbo suppressed 
processes. 



various J /if) ^ D*J^ + M modes are presented and comparisons of our numerical results with 
that estimated in other theoretical models are also investigated at length in this section. The 
final section is devoted to the discussions and conclusions. It is noted that since most of the 
form factors applied in this work were obtained in our previous work, we generally refer the 
readers to it for some details of the derivation and how to achieve the numerical values. 



II. NON-LEPTONIC DECAYS J/^j D^,) + M IN FACTORIZATION APPROACH 

For the non-leptonic weak decays of J/ip —>■ D(^s) + M, the standard method is integrating 
out the heavy PI^— boson and obtaining a low energy effective Hamiltonian for c quark decay 
which is given by 

Heffic ^ quq') = (CiQi + C2Q2) , (1) 
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where g(g') represents the down type quarks s and d; Kg(Kjg') are CKM matrix elements; 
and the operators Qi, Q2 are respectively 



Ql = Qalt^il - 75)Ca g^7^(l - 75)M/3, Q2 = Qa^f^il " 75)C/3 QpYi^ " l5)Ua- (2) 

It should be pointed out that the penguin operators are neglected in this work due to the 
smallness of Wilson coefficients for such operators, which also indicates that CP symmetry 
is well respected within the accepted assumption. 

With the free quark decay amphtude, we can proceed to calculate the transition ampli- 
tudes for J/ip D+M at hadron level, which can be obtained by sandwiching the free-quark 
operators between the initial and final mesonic states. Consequently, the hadronic matrix 
elements {DM\Qi\J /ip) which depend on the strong interactions need to be computed. The 
evaluation is indeed the main challenge in the heavy flavor physics due to our poor knowl- 
edge with respect to the non-perturbative QCD. Owing to the painstaking efforts in theory, 
several systemic approaches for hadronic B decays have been explored based on the expan- 



sion in small parameters 20]. However, a systematic theoretical method concerning the 
open-charm decays is still not available yet due to the fact that the accessible charm quark 
mass is not so heavy in reality. As the ffist order approximation, we may be able to apply 
the vacuum saturation approximation to factorize the four-quark operator matrix elements 
l^DM\Qi\J/il)). The consistency of the theoretical prediction with data (may be available 
in the future) will serve as an examination of such approximation in the heavy- quarkonium 
system as mentioned in the introduction. To be more specific, the factorization ansatz [14j 
states that the matrix elements can be factorized into a product of two single matrix ele- 
ments of currents (M| Ji|0)(D| J2|^/^) where one is parameterized by the decay constant of 
the emitted light meson and the other is represented by the form factors responsible for the 
transition of J/ip into the recoiled charmed meson. 

The decay constants for pseudoscalar (P) and vector {V) mesons are defined as follows 

(P(g)|A^|0) = -ifpq^, 
iy{q,e)\V,\Q) = fvmve;, (3) 

where the axial vector current represents gi7^75g2 and the vector current represents 
'li1fj.Q2', e is the polarization vector of V. The matrix elements (D|g7^(l — '-)^)c\J /ip) are 



parameterized in terms of various form factors as 
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where q = pi—p2 and the convention T^T[^^'~iulpl(jlb] = 'H^^iupa is adopted. For the transition 
of J/-?/^ into a charmed pseudoscalar meson which is induced by the weak current, there are 
four independent form factors: V, Aq, Ai, A2] while there are ten form factors for J/ip 
transiting into a charmed vector meson which are parameterized as Ai{i = 1, 2, 3, 4), Vj{j = 
1,2,3,4,5,6). 

According to the quark diagrams in Fig. [H the decays are classified into two categories: 
color allowed and suppressed processes. For the color allowed processes, the decay amplitudes 
are proportional to 



(6) 



with Nc being the color number of QCD. Because Ci ~ 1 and C2 ~ a^, ai is estimated to 
be of order 1. While for the decays shown in Fig. [T]^b), the amplitude is proportional to 



a2 = C2 + Ci/iV,. 



(7) 



As 02/01 ~ i/Nc, this kind of decays are usually named as color-suppressed processes. 

For the decays of c ^ sud which is the Cabibbo-favored process, the CKM element VcsVud 
responsible for these modes is close to 1. For the Cabibbo-suppressed transitions of c — dud 
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and c sus, the corresponding CKM parameters VcdVud and VcsVus are suppressed by a 
factor sin6'c ~ 0.22 with 6c being the Cabibbo angle. The doubly suppressed processes, such 
as c — i> dus which are suppressed by siv?6c, are neglected in our case. Thus, the prevailing 
decay modes are both color allowed and Cabibbo favored ones. The less dominant modes are 
the color suppressed but Cabibbo favored or color allowed but Cabibbo suppressed processes. 
These are the processes we will focus on in this study. Note that there is no annihilation 
type contributions in our case at all. 

Now, we are able to write down the decay amplitudes associating with the non-leptonic 
two-body decays of J/ip explicitly based on the information we achieved before. In light of 
the characters of final states, three different types of processes J/ip DP, DV , D*P and 
D*V will be investigated one by one in the following sections. As for J/tp decaying into two 
pseudoscalars where one is a D meson and the other is a light meson P, the decay amplitude 
is written as 



A{j/ij ^ DP) = {DP\neff\J/ij) 



■■V*Vuq'ai 2m^(e^ ■ q)fpAo{q^), 



s// 1 " / r / — ^2 ' eg ■ 

where q denotes the momentum of light emitted meson; a, is the effective coefficients with 
tti for color allowed process and 02 for color suppressed process. 
The decay amplitude of J/ip ^ DV decay is given as 



G 



A{J/^p DV) = ^V*gVuq'a^ fvmv \ -et,uape*y^e'^P^p'^D—^ 
V Z [ Tn-ii) + 



2V{q^ 



+i{m^ + mD){e^ ■ ey)Ai{q ) + — - — — 2^2(5 ) 



+ rriD 



A{J/i;-.D*P) = t^V:^V^q,a,fp\2te,,^pp'tp^,e^etA,iq') 



(9) 



the decay amplitude of J/ip D*P decay can be shown as 

. G p 
72''^ 



[ml-ml,)V^{q')-q'V2{q' 



- nil') + V^iq') - v,{q' 



(10) 



Lastly, the expressions for J/ip D*V can be readily derived from Eqs. ([3],|4]) as 
Gp 



A{J/i) D*V) 



72 



,0 *I3 
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(e^-eD0[4-bi+P2)V^i(g')] + 



- (e^ ■ q)e*^, ■ e*yV,{q') + (e^,, ■ q)e^ ■ e^yVeiq')^ (11) 
where ey denotes the polarization vector of hght emitted mesons. 

III. DECAY RATES FOR NON-LEPTONIC WEAK DECAYS OF J/tp 
A. Input parameters 

The decay rates of the non-leptonic decays J/ifj ^ D + M are written as 



T^^nM = \^ \A{J/^ ^ DM)f (12) 



where pd denotes the three-momentum of the final D meson in the rest frame of J/tp and 
the factor "|" is due to the spin average of J/ip. In order to calculate the decay rates, 
the input parameters including the CKM parameters, effective Wilson coefficients, decay 
constants and transition form factors are necessary. The CKM parameters are taken from 



ref. 



2l| 



Kd = 0.974, Ks = 0.227, = 0.227, = 0.973. (13) 

The effective Wilson coefficients are determined as 

ai = 1.26, a2 = -0.51, (14) 

which are extracted from the isospin analysis for D Kit decays with the help of the 
factorization ansatz [22]. 



The decay constants for light mesons are taken as 



21 



23| 



= 0.131 GeV, /k = 0.160 GeV, 

fp = 0.209 ± 0.002 GeV, fx* = 0.217 ± 0.005 GeV, (15) 

where the pseudoscalar decay constants are determined from the combined rate for P l^ui 
and P l^vn experimentally and vector meson longitudinal decay constants are extracted 



from the data on r ^ (p , A"* )i/t- |21 |. 
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TABLE 
model 


] I: The form factors Aq and Ai at 
3] and QCDSR approach. 


= responsible 


for the decays oi J/ip - 




Models ^^^^ 


^0 






BSW 


0.61 


0.66 


0.68 


0.78 



QCDSR 0.27 0.37 0.27 0.38 



TABLE IL Branching ratios of non-leptonic decays of J ftp — > (in units of 10 ^^). 







other works 


this study 


BR{J/^ - 




17.4 [4j 
10.0 [5] 


2 0+°-^ 


BR{J/%b - 




1.10 [4] 


16+°°2 

'-'•-'^"-0.02 


BR{J/iP - 


Dn) 


1.10 [4] 


o.o8ol[J:°2 


BR{J/^ - 


DK) 




36+°-^° 



Besides, the values of all the revelent transition form factors Fi{q^) are taken from our 
earlier study 6| where the detailed expression are presented. In the literature, the form 
factors Aq and Ai have been calculated by various authors {4, 5|, which are grouped in the 

We can see that the 



Table [H together with the numbers obtained in the QCD sum rules 
form factors at zero momentum transfer predicted in the BSW model 
greater than that in the QCD sum rules by a factor 2. 



24| are approximately 



B. Branching ratios of non-leptonic decays 

The numerical results of branching ratios for non-leptonic decays of J/ip D{s)P are 
presented in Tables [Tll where the numbers obtained in ref. 0, are also collected together 
for a comparison. Here, the results are given for decays which including the charge conjugate 
process, for instance, BR^J/ip —>■ Dgir) is the branching ratio for decays oi J/ip Dfir^ + 

Table HI] shows that the decay rate for color allowed and Cabibbo favored channel J/ih 

n 

DsTi calculated in this work is five times smaller than that given in Ref.[5|. Such discrepancy 
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may be attributed to two aspects: Firstly, an SU{4)fijk rotation matrix is employed in 
Ref. to relate the J/ip ^ Ds transition to D K* decay, and the form factor Ao{0) 
is estimated as 0.7 ~ 0.8, which is almost twice as that computed in the QCD sum rules. 



Secondly, the experimental data on total decay width o: 
however, this value has been updated to 93.4 ±2.1 keV 



2l|. 



J/x/j used in Ref. js] is 67.0 keV 



As for the Cabibbo suppressed but color-allowed mode J/ip ^ DgK, the following relation 



Ri 



BR{J/iP D,K) 
BR{J/%1) ^ D.-k) 



Vusf. 



K 



VudflT 



0.081 



(16) 



is achieved in the factorization assumption. Similarly, we can define a parameter R2 as 



R2 = 



BR{J/i) Dtt) 



0.032, 



(17) 



BR{J/ij DsTt) 

which is also in agreement with that listed in Table [Tll as long as the phase space is properly 
considered for these two channels. 

Now, we move on to the discussions of color suppressed mode J/ip ^ DK. A ratio of 
decay rates between it and Jjil) ^ Dgir can be estimated as 



R^ = 



BR{J/iP DK) 



0-2 



An 



m 



K) 



0.18, 



BR{J/iP Dsn) 

which is consistent with that collected in Table HTl In addition, we should emphasize that the 
ratio i?3 is quite sensitive to the effective Wilson coefficient 02, which can receive considerable 

3, 



corrections 



251, 



261, 



23, 



29 



30 



31 



32l | due to uncertainties of the renormalization scale. 



higher order effects together with non-factorizable contributions, where we also refer to [20] 
for a recent comment. 

Table UTTl collects the numerical results for J/ip ^ ^(s)^ -^(s)^ decays. As one can 
see, the branching ratio of J/t/j D^p computed in j^j] is 5.8 times larger than that evaluated 
in this work. It is shown in [4], the dominant contributions for decay width oi B ^ D^*^P 
are from the form factor Ai{q^) corresponding to the S partial wave in the final states. As 
listed in Table. [H the number of form factor Ai derived in the BSW model is 2.1 times 
greater than that in terms of the QCD sum rules, which can indeed result in an enormous 
discrepancy for the branching fraction oi J/ip —>■ Dg obtained in two different approaches. 
Moreover, the ratio 

BR{J/^ ^ Dsp) 



Ra 



(19) 
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TABLE III: Branching ratios of non-leptonic decays of J ftp ^{s)^ ^^'^ ^(s)^ units of 10 ^^). 







other works 


this study 


RR( Jill) - 


n n\ 




in f,+3.0 






4 94 f4l 


n S9+0-22 




r)n\ 


4 4n \d\ 


n 49+0-18 
•-•■^^-o.os 


BR{J/iP - 


^ DK*) 




1.54+[!i 

— U.oo 


BR{J/iP - 


- D» 


— 


15.0lJi 


BR{J/iP - 


- DtK) 




1 1+0.08 
—0.04 


BR{J/iP - 


D*tt) 






BR{J/iP - 


^ D*K) 




2.61°:^ 


TABLE IV: Branching 


; ratios of non-leptonic decays of J/V' - 


L>; ,y (in units of 10"^°). 


Channels 




this study 


BR{J/ij - 


- Dip) 






BR{J/iP ^ 


> DIK*) 




^•D_o.4 


BR{J/^ - 


-D*p) 




9 O+0.6 


BR{J/tl; ^ 


. D*K*) 







is usually introduced from a viewpoint of experiment, whose value is estimated as 6.3 and 
4.2 respectively in the framework of QCDSR and BSW model. Therefore, the decay of 
J/ip DsP is more detectable than the corresponding pseudoscalar channel J/ip ^ Ds'k in 
experiment. Moreover, it can be seen that the ratio of decay rates is not sensitive to the 
absolute magnitude of the transition form factors on account of the large cancelations of the 
non-perturbative effects. The relative magnitude of decay widths for the Cabibbo-suppressed 
as well as color suppressed processes to the mode of J/ip ^ D^p can be readily derived by 
following the discussions on J/ip ^ ^{s)P and will not be repeated again. 

Furthermore, we group the decay rates for dominant channels of J/ip —>■ D*^-^V in Ta- 
ble IIVI from which we can observe that BR{J/ip D*p) is as large as 5.3 x 10~^ and stands 
as the most promising mode to be measured at BESIII. Such finding presents a striking con- 
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trast to the argument given by the authors of Ref. [l| where the authors claimed that a 
specific non-leptonic decay channel like J/if) ^ D^*'>M (M = vr,p...) is hardly to be detected 
owing to the tiny branching factions for these processes. In addition, it is also helpful to 
define the following two ratios 

^ BR(J/^ ^ ^ BR(J/^ ^ D^p) 

' " BR( J/^ D,7r) ' ' " BR( J/^ ^ D,p) ' ^ ^ 

which characterize the relative size of branching fractions to distinguish the final states with 
vector and pseudoscalar ones respectively in the non-leptonic two-body weak decays of J /ip. 
The numbers of and Rq are evaluated as 7.5 and 4.2 in the QCD sum rules, while they 
are determined as 3.5 and 1.4 respectively with the ISGW model in the framework of heavy 
quark spin symmetry [Ij. Such discrepancies can be attributed to the different values of form 
factors employed in the numerical calculations. 
Moreover, we also mention that 

^ BR(J/V; ^ D:K*) BR(J/^ Dp) 
^ BR{J/tl) ^ D*p) BR{J/iP ^ DsK*) ^ ' 

should be equal to 1 in the heavy quark limit. However, this ratio is estimated as 0.48 in 
the QCD sum rules owing to a serious suppression factor from the phase space for the decay 
of J/ip D*K* for the limited charm quark mass. 

Combining the Table [Tll IIIII and IIVI we find that the branching ratio for inclusive weak 
decay of J/ip can be as large as 1.3 x 10^^, which is also in remarkable agreement with the 
naive estimation 

BR{J/i) ^ + ...) ^ ^^-^ ^ 1.4 X 10~l (22) 

IV. DISCUSSIONS AND CONCLUSIONS 

Since J/ip mainly decays via strong and electromagnetic interactions, its weak decays 
usually take small fractions which cannot be measured by available experimental apparatus. 
On other aspect, however, because J /if) contains two heavy constituents, its weak decay may 
possess a unique character. Indeed weak decays of J / ip may offer an ideal platform to examine 
the mechanism which governs the hadronization process, without possible contamination 
from the light spectator as well as one may determine its fundamental parameters such as 

12 



the CKM matrix which can be a complementary test to the values obtained in D decays. 
It is lucky for high energy physicists that a tremendous database on J/ip will be available 
in the forthcoming BESIII and the measurements on the weak decays of J/iJj may become 
possible. 

As is well known, the essential challenge in the theoretical calculations on the rates of 
weak decays of J/ip is to disentangle the underlying weak- interaction transitions from the 
notorious effects owing to strong interactions reasonably [343 • previous paper (gI , the 

transition form factors in the semi-leptonic weak decays of J/ip have been investigated to 
the leading order of as based on QCD sum rules, where the non-perturbative QCD dynamics 
is characterized by a few universal parameters. The branching ratios for dominant exclusive 
processes are evaluated and their order of magnitude is typically at 10~^°. Obviously based 
on the factorization assumption, the form factors obtained for the semi-leptonic decays can 
be applied to study the non-leptonic decays. 

This paper can be viewed as a continuation of our earlier work j^. We present a com- 
prehensive study of non-leptonic decays of J/ip —>■ D(^s) + M based on the factorization 
assumption and apply the transition form factors calculated in the QCD sum rules. It is 
observed that the sum of the branching fractions for the dominant non-leptonic decays of 
J/ip ^ D^TT, D~p, D*~iT, and D*~p as well as their charge conjugate channels can reach as 
large as 0.82 x 10~®, a special decay mode J/ip ^ D*p can even arrive at 5.3 x 10~^, which is 
hopefully to be marginally detected in the e~^e~ colliders in view of the large database of the 
BESIII. Our results are in agreement with the finding in Ref.|l| that J/ip decays to vector 
charmed meson D*~ more favorably than to the pseudoscalar one, however, the ratios of 
these two channels calculated in this work is twice or three times larger than that given by 
Ref. where heavy quark spin symmetry and the non- recoil approximation were adopted 
and the ISGW model was employed to compute the single form factor rii2. 

As there is a light to see a possibility of measuring weak decays of J/ ip which have obvious 
advantages for getting insight to the physics picture, we strongly urge our experimental 
colleagues to search for vector charmed mesons productions in J/ip decays at BESIII 



Moreover, from the theoretical side, it shou 
tions for the heavy quarkonium system 



35 
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De_emphasized that Coulomb-type correc- 
are not included in the computations 



of form factors in the QCD sum rues, which could induce additional uncertainties to the 
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evaluation of the branching fractions for non-leptonic two body decays of J/ip. However, 
one can still trust the order of magnitude gained in this work, since while calculating the 
form factors which need to deal with the three-point correlations, most uncertainties origi- 
nating from Coulomb-like corrections are canceled by that in the two-point correlations for 
evaluating the decay constant of J/ip. Apart from weak decays of J/ip presented in this 
paper, weak decays of T are suggested to be explored seriously in a complementary fashion 
from both the theoretical and experimental point of view. 
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